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ABSTRACT 

We present the results of two high-statistics studies of the gluon propagator in 
the Landau gauge, at /5 = 6.0, on different lattice volumes. The dependence 
of the propagator on the momenta is well described by the expression G{k^^ = 

\M^ + Z-k^{k^ I K^y^)]^ \ We obtain a precise determination of 7/ = 0.532(12), 

\l and verify that M"^ does not vanish in the infinite volume limit. 

On 

a^ 

^ The non-perturbative investigation of the behavior of the basic fields of the QCD 

Q Lagrangian is crucial to shed light on the mechanism of confinement and can be 

Q achieved through numerical lattice computations of the Gluon Propagator-^ [GP]- 

The Euclidean GP in the Landau Gauge is: 

<N 



D,,{k) = j d'xTT{A,{x)A,{Q))e-'^- = G{k') ( ^ " ^) , (1) 



where /x, = 1, . . . , 4 and the trace is intended over color indices. 

Recently, there has been much effort in trying to obtain a non perturbative form for 
G{k^), both from analytic'^"'' and numerical^" -^-^ analyses. With the present studies, 
we investigate the non-perturbative form of G(k^) and its behaviour in the infinite 
^ volume limit. 

In tab.l, the parameters of our simulations are summarized. The gauge fields 
have been generated with a Hybrid Monte Carlo algorithm^^^. The Landau gauge- 
^ fixing has been performed and checked carefully, being this a crucial point when 

4^ dealing with gauge-dependent quantities. The fluctuation left-over after gauge-flxing 

{d^A^(x))i,^tt < 10~^, are absolutely negligible-^, with respect to the statistical errors. 

This can be checked since the condition d^A^(x) = implies dtAo(0,t) = 0. One can 
deflne the correlation {Ao(t)Ao(0)) at zero momentum, and study its time derivative. 

We have shown-^'^ that it is zero within errors, and the gluon fleld Ao(0, t) is constant at 
the level of 0.008% on each individual conflguration. These results demonstrate that 
the present gauge-flxing procedure is the most effective, among those implemented in 
the literature^" . 

Accordingly to eqn.l, we compute 2-pt functions of the gluon fleld, deflned in 
term of the link variable Ufj_{x) as Afj_{x) = [Ufj_{x) — u}i,{x)]/2i, which, using spectral 
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/5 


^ confs. 


Volume 




6.0 


1000 


16^ X 32 


< 10"^ 


6.0 


500 


24^ X 48 


< 10"^ 



Table 1: Summary of the parameters of our simulations. 



decomposition and translation invariance, can be written as 



D{t,k) = Y.i:T{A,{x)A,{Q))^ 



-E,t 



J = 1, ... ,3, 



(2) 



where the sum is over the states which couple to the gluon field and Ei is the energy 
of the state From eq. (2), the effective energy, defined as 



uJesit, k) = log 



D{t,k) 
D{t + l,ky 



(3) 



should be a decreasing function of the time, for any value of the momentum k. In fig. 
1, we show that a;eff(t, k) is increasing with time, for all the momentum combinations 
considered. Hence, it is impossible-^'^ to fit the GP to a sum of single particle pole 
function, neither if physical states (J\fi > 0) nor "ghost" (J\fi < 0) are considered. 
This is an unacceptable feature for the propagator of a physical particle. To avoid 
systematic uncertainties due to the presence of infrared and ultraviolet cut-offs on 
the lattice, one can study the GP in the intermediate region of momenta. We find 
that G(k^) is well described by the following modeling function 



G{k') 



1 



on the lattice 



M2 + 



(4) 



which depends on the parameters r/, and Zl. The stability and the quality of the 
fits have been checked in different ways-^'^. On the two volumes, we find (see fig. 2) 



V = W 



32 



MI-- 
Zl = 
V = 

Xndof 



= 2.8(1) : 
9.01(4) 

0.56(6) 
= 1.5 



10" 
: 10" 



F = 24^ X 48 <^ 



V = 
Ml 
Zl- 

V = 

Xndof 



24^ X 48 
= 4.46(9) : 
= 0.102(1) 
0.532(12) 
1.08 



10" 



(5) 

We observe that the finite volume has a very small effect in the value of the 
anomalous dimension, provided that the range of k^ is large enough. Indeed, we obtain 
a fairly accurate determination of rj. The two determinations of are inconsistent, 
and increases with the volume. This feature rules out the hypothesis that the 
non-zero value of is merely due to finite volume effects. If this were the case, we 
would expect to scale roughly as 1/ . It is possible to try a first, very crude 
extrapolation of the value of to the infinite volume limit. Using 



M\V) = M\V = oo) + cost 



we get M\V = oo) = 6.202(8) x 10"^ (6) 




Figure 1: Effective energy for gluon 2-pt functions. The four curves correspond to 
the following momenta: □ : A; = 0, O : A; = (27r/24, 0, 0), x: k = (27r/24, 27r/24, 0) and 
: k = (47r/24, 0, 0). Figure 2: Best fit of the propagator in momentum space to the 
function 5. The figure corresponds to the case V = 24'^ X 48. 



Using a ~ 2GeV, as determined by several simulations at /5 = 6.0, we obtain 

Mp\^3 ~ (160 MeVf ~ {AqcdY. (7) 

An interpretation of this result, in connection with colour confinement is, at 
present, absent. However, we stress that it is fundamental to understand the be- 
haviour of both 7/ and in the continuum limit. 
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